finish
/clear
/filn, E322  
/title, FEA of connecting zone of nozzle to cylinder
! *************趨*******************
Rci=1000	                !Ͳڰ뾶
tc=30	                        !Ͳ
Rco=Rci+tc                      !Ͳ뾶
Lc=4000	                        !Ͳ峤
Rno=530	                        !ӹ뾶
tn=15                           !ӹܺ
Rni=Rno-tn	                !ӹڰ뾶
Hmin=sqrt(Rci**2-Rni**2)        !ӹС쳤
Lni=50			        !ӹ쳤
Ln=500	                        !ӹ쳤
rr1=30                          !Բǰ뾶
rr2=15	                        !ڲԲǰ뾶
pi=1.2	                        !ѹ
pc=pi*Rci**2/(Rco**2-Rci**2)	!Ͳ˲ƽغ
pn=pi*Rni**2/(Rno**2-Rni**2)    !ӹܶ˲ƽغ
!****************ǰ***************************
/prep7
et,1,45	                        !嵥Ԫ
mp,ex,1,2e5	                !ϵĵģ
mp,nuxy,1,0.3	                !ϵĲɱ
!****************ģ***************************
cylind,Rco,Rci,0,-Lc/2,90,270,	!Ͳ
wpoff,0,0,-Lc/2	                !-ZƶLc/2
wprot,0,90, 	                !yzת90o
cylind,Rno,Rni,-Ln-Rci-tc,-Hmin+Lni,90,180,      !ɽӹ
vovlap,all	                !overlap
vdel,3,,,1              	!ɾͲϿʧȥĲ
vsel,s,,,7              	!ѡͲ 
*afun,deg	                !趨ǶȺеλΪǶ
ang1=2*nint(asin(Rno/Rci))	!ӹиǶ
wprot,0,0,-90+ang1	        !תϵ
vsbw,all                	!иͲ
afillt,21,12,rr1          	!ͲӹԲ
al,33,60,5	                !Բ
al,37,61,7	                !Բ
va,12,14,16,34,35	        !Բ
afillt,23,13,rr2	        !ͲڱӹԲ
al,24,40,68	                !Բ
al,25,41,69	                !Բ
va,13,23,31,39,40	        !Բ
alls,all	                !ȫѡ
askin,62,70	                !ԲǱ߽Ƥи
vsba,2,36	                !иͲ
wprot,0,,90-ang1        	!תϵ
wprot,0,90,0	                !תϵ
wpoff,0,0,-2*Rni        	!ƶϵӹиλ
vsbw,9	                        !иͲ
vsbw,1                      	!иͲ
!*******************************************
lsel,s,loc,y,Ln+rci+tc	        !ѡӹܶ˲߶
lsel,u,length,,0,tc	        !ȥΪ0tc߶
lesize,all,,,20,,1           	!趨ѡ߶εʷ
alls,all
lesize,60,,,4,,1        	!趨Բĳһ߶εʷ
lesize,45,,,12,,1               !趨ڽӹֶܷ
lesize,57,,,12,,1
asel,s,,,45,46,1,1	        !ѡͲϽӹӰ
accat,all	                !
alls,all                        !ȫѡ
vsweep,all                      !ӳ仮
alls
nummrg,all              	!ϲͬ
numcmp,all	                !ѹ
fini
!*******************************************
/solu
asel,s,loc,z,0	                !ѡͲ
SFA,all,1,PRES,-pc	        !ʩӶƽغ
asel,s,loc,x,0 	                !ѡԳ
asel,a,loc,z,-Lc/2	        !ѡԳ
DA,all,SYMM	                !ʩӶԳԼ
asel,s,loc,y,Ln+Rci+Tc	        !ѡӹܶ
DA,all,UY, 	                !Լλ
asel,s,,,6	                !ѡڱ
asel,a,,,7
asel,a,,,15
asel,a,,,17
asel,a,,,19
asel,a,,,23
asel,a,,,29
asel,a,,,42
asel,a,,,47
asel,a,,,50
SFA,all,1,PRES,pi	        !ʩѹ
alls
solve	                        !
fini
!*******************************************
/post1
PLNSOL,S,INT,0,1		!ʾӦͼ
fini
